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ABSTRACT 

hi Part I of this paper, Popov's Theorem PI is first introduced and then 
a new Theorem I is formulated. The proof is given in Appendix I. An illustrative 
example shows that the result obtained from Theorem I agrees with that obtained 
from Lur'e’s Theorem. In Part n the linear part transfer function may have 
poles along the imaginary axis with real positive residues. The nonlinear func- 
tion f(e) is bounded as well as continuous. Popov's Theorem PH is extended to 
form a new Theorem n, which gives the condition for quasi-asymptotic stability. 
Two corollaries are also given. Corollary Ha gives the condition for asymptotic 
stability. The proof of Theorem n and its coroUaries is given in Appendix EL 
Three examples check with the results of analog computer studies. 
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FOREWORD 


On April 19, 1966, at the Space Science Seminar at NASA's George C. 
Marshall Space Flight Center, Dr. Y. H. Ku of the University of Pennsylvania 
presented a lecture entitled Nonlinear Oscillations in High-Order Systems. This 
presentation was based largely on a paper published in the Journal of The 
Franklin Institute, June, 1965. 

This Technical Memorandum is a reprint of the original Franklin Institute 
paper, Extensions of Popov's Theorems for Stability of Nonlinear Control Sys- 
tems. The seminar report also included 30 or 40 additional figures of cyclic 
stable and unstable functions. 
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EXTENSION OF POPOV'S THEOREMS FOR STABILITY 
OF NONLINEAR CONTROL SYSTEMS 

INTRODUCTION 


Popov has recently given two theorems by which we can determine the 
stability of nonlinear control systems. Theorem PI gives the asymptotic 
stability of the system involving a nonlinear function f(e ) by the condition 
Re > 0, where G(jw) denotes the linear part transfer func- 

tion, and q is a non-negative constant. The new Theorem I is an extension of 
Theorem PI such that the restriction 7 > 0 is now changed to 7 > 0. Similar 
to Theorem PII, the new Theorem II is applicable to control systems with its 
linear part transfer function having pairs of imaginary poles. While Theorem 
PII gives the condition for hyperstability, Theorem II gives the condition for 
quasi-asymptotie stability. Corollary Ila then gives the condition for asymp- 
totic stability whenever q can be set to zero. 


EXTENSION OF POPOV'S FIRST THEOREM 


Popov's Theorem PI (1) is first introduced in this section. Extension of this 
Theorem will be given as new Theorem I. (The proof of Theorem I is given in 
Appendix I.) An example shows that the result obtained from Theorem I 
agrees with that obtained from Lur’e's Theorem. 

Popov's Theorem PI: If a non-negative quantity q exists such that for 
all real w, the inequality 

ReC(l+i^)(?(ico)]>0 (1) 

takes place, then the trivial solution of the system is asymptotically stable 
provided a given set of assumptions are valid. 

Kalman (2) restated Popov's Theorem PI as follows: 

"Assume that A is stable and that 7 > 0. Then the system is g.a.s. if 
the condition 

Re { {2ay + jo>q)lC'(jwI - A)" l b + 7 //«]} > 0 (2) 

for all real « holds for 2a7 = 1 and some q > 0." 

In the above quotation, we have replaced the original symbol F by A for a 
real n X n matrix, h by C for a real n-vector, g by b for another real n- vector, 
and the constants p and 0 by 7 and g, respectively. The prime denotes the 
transpose. Comparison of the two expressions gives the transfer function of the 
system linear part as 



G(joj) = C'(jwl - A)~*b + y/ja> (7 > 0). (3) 

Note that on the righthand side of Eq. 3, the first term contains the poles with 
negative real parts while the second term denotes a pole at the origin. 

The new Theorem I tries to remove the restriction that 7 > 0 so as to 
include the case 7 = 0 in its application to feedback control systems shown in 
Fig. 1 and introduces a new constant h which is non-negative. The system can 
be defined by the following differential equations : 

x = Ax + bf(e), (4) 

e = — C'x — hf(e ), (5) 

where x denotes a real n-vector, x = dx/dt , A is a real n X n matrix whose 
eigenvalues may include zero, b and C arc real n-vectors, and h is a non- 
negative constant. The nonlinear function /(e) is a continuous function of e f 
with the conditions : 

0 < ef(e) < ke*(e ^ 0); /( 0) - 0. (6) 



Dls) * $ 3 + IOs 2 + 31s ♦ 30 
N(s)» E + l 


Fig. 1. Block diagram of a nonlinear Fro. 2, Block diagram of example 1. 

control eyBtem. 

It is further assumed that f'(e) = df(e)/de is continuous, and 1 + hf f {e) ^ 0. 
Referring to Fig. 1, the linear part transfer function is defined by 

G 1 (s ) = C'(sl - A)-*b + h. (7) 

For s = jo) y Eq. 7 gives Gi(jw), where 7 is absorbed in the first term on the 
righthand side. The new Theorem I can be stated as follows: 

Theorem I: If a non-negative quantity q exists such that for all real w 
the inequality 1 

l/k + Re [(1 + juq)Gi(jw)~\ > 0 (8) 

takes place, then the trivial solution of the system Eqs. 4-5 is asymptotically 
stable provided the assumptions made are valid. 

Substituting Gi(ja)) into Eq. 8 gives 

l/k + Re [(1 + i«$)C'(j«I - A) -I b] + h > 0. (9) 


1 Note the factor (1 /k) appeared in Popov and ITalanay (6) and in Aizerman and Gant- 
macher (7), p. 52). However the expression Eq. 8 in the present paper refers to G\(jai) which 
is more general than G(ja>) in (6), or W (jot) in (7). Popov (4) did introduce the factor h (but 
not 1/jfe). Maigarin (8) gave in (2.14) an additional term N which corresponds to h in Eq. 7, 
but the (1 /k) factor wasjnot included. 
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When h * 0, the condition that /'(e) is continuous can be Temoved and that 
the condition 1 + hf 0 is automatically satisfied. The proof is given in 
Appendix I. 


Example 1: Given a control system as shown in Fig. 1 with the linear part 
transfer function: 

GlW = (s + 2)(* + 3)(# + 5) ’ fc = 0; 7==0 - (10) 

The left hand side of inequality E<p 8 gives 

1/Jr + Re [(1 + jo9q)Gi(j<a)] = + A* o 4 + Ai<a* + Ao)/B(od), (11) 

where A t =l/k, A t =9q- 1+38/Jr, A^g+21+361/*, A 0 =30(l+30/*), and 
B(<a) = 100(w* — 3)* + ^(w 1 — 31)*. Expression Eq. 11 is greater than zero 
when q > 1/9 for any positive * and any «. Hence the system is asymptotically 
stable , according to Theorem I. 

This result checks with the result by LurVs Theorem as shown in (3). 


For this example, a set of state variables x can be chosen as xi, x*, and x« 
where X\ is shown in Fig. 2. Referring to Eqs. 4 and 5, we get 



0 1 0 


0 


t 

A - 

0 0 1 

b - 

0 

c - 

i 


-30 -31 -10 


1 


0 


For h = 0, substituting Eqs. 12 into Eq. 7 gives (?i(«) specified by Eq. 10. 
Note that from Fig. 2, we get the output c, which is the negative of e in an 
autonomous system (with input r = 0), as 

c ■» Xi + ±\ — X\ + Xj. (13) 

From Fig. 2, we also get the differential equation 

f(e) « x, + lOx, + 3 lx* + 30x, (14) 

where Xi = x a and Xi = Xi are defined in matrix A. The variables x i, x t , and x$ 
are known as phase-space variables. 


EXTENSION OF POPOV'S SECOND THEOREM 


Popov's Theorem PII (4) is first introduced in this section. Extension of 
this Theorem will be given as new Theorem II. (The proof of Theorem II is 
given in Appendix II.) Three examples are given to illustrate the application 
of Theorem II. Analog computer results are presented to check the stability 
of the three control systems studied. 


Popov 9 s Theorem PII: 

In order for the solution x — 0 of the non-degenerate system Eqs. 15-16 
to be hyperstable it is necessary and sufficient that the transfer function C?i(a), 
defined in Eq. 18, be a real positive function. 

The system considered in Theorem PII can be defined by the equations 

i « Ax + b*(<) (15) 

0(0 - C'x + Hit) (16) 

which satisfy the following inequality for any T > 0 : 

f 0(00(0* < 8 sup |x(0|, (17) 

JO 0<t<T 


where & is a non-negative constant. The transfer function <?i(s) is defined by 


GiM = i + g«m. as) 

where Gi(s) is defined in Eq. 7, p< > 0 and co< are constants. The solution of 
the system Eqs. 15-16 is hyperstable if there exists a constant K such that any 
solution of the system satisfies the inequality 

|x(0 < K(|x(0)| +5) (19) 


for all 0(0 that satisfies inequality 17. 

For the sake pf developing new Theorem II, let tf>(0 = f[e{t)~\ and 0(0 
= — e(t) with the conditions given in Eqs. 6. Substituting these into Eq. 17 
gives 




ef{e)dt<h sup |x(0|. 

o<t<r 


( 20 ) 


This is satisfied for any non-negative 5 (including 6 « 0). The system Eqs. 
15-16 now take the form of Eqs. 4-5, where the A matrix has eigenvalues on 
the imaginary axis as well as on the left half of the s-plane. 

The new Theorem II Can be stated as follows: 


Theorem II: If /(e) is bounded and a non-negative quantity q exists such 
that for all real o) the inequality 

l/k + Re [(1 + i«7)G 0 (i«)] > 0 (21) 

takes place, then the trivial solution of the extended system Eqs. 4-5, with 
its transfer function (?i(s) defined by Eq. 18, is quasi-asymptotically stable 
provided the above assumptions are valid. 

Note that the definition of quasi-asymptolic stability follows that given in 
(5). The proof is given in Appendix II. 

Corollary Ila: If q can be set to zero, then the trivial solution of the 
extended system Eqs. 4-5 is asymptotically stable provided the assumptions 
are valid. 
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Corollary lib: If q can be set to zero, /(e) is continuous, but the restric- 
tions on /'(e) are disregarded, then the trivial solution of the extended system 
Eqs. 4-5 is stable. 

These corollaries are discussed in Appendix II. 


Example 2: Given a control system shown in Fig. 1 with the linear part 
transfer function : 


Gi(s) = 


3s 4 + 4s* + 5s 5 + 8 + 1 
s(3s + l)(s* + 1) 


( 22 ) 


and /(e) = tanhe, which is bounded and continuous. The derivative /'(e) 
= seen 5 e is continuous. Note that/(0) = 0 and i + kf'(e) = l + sech* e > 1 
which is not equal to zero. Dividing the numerator of Eq. 22 by its denominator 
gives 5=1 and 




3s* 4- 2s* + 1 
s(3s + l)(s* + 1) 


(23) 


Expansion into partial fractions gives 


Gi («) 


s i 3 
8* + 1 + 8 38 +1 


+ 1 . 


(24) 



Comparison with Eq. 18 shows that pi = 1 and = 1. Equation 24 also 
shows that the coefficient of 1/s is y = 1, and this term denotes a pole at the 
origin. Note that for this example, (?o(s) in Eq. 18 is given by 



3 . . _ 3s* + s + l 

3s + 1 + 1 _ s (3s + 1) • 


(25) 


Referring to inequality 6, we get 

0 < e tanh e < fre*(e 7 ^ 0) ; tanh 0 = 0 (26) 


and k can be chosen as 1. Then applying inequality 21 of Theorem II gives 


1 + Re [(1 +M) + l)] = (Q I+lL— - ° 


for 


q > 1. (27) 




According to Theorem II, the system is quasi-asymptotically stable. The 
system can be represented by the block diagram shown in Fig. 3. The variables 
chosen for analog computer study are: x 0 , x it zj, and z*. The system 
equations are : 


and 


i« = /(c) 

*1 = — ixi + /(e) 

Zj = z, 

i\ “ - Wl S ?l + /(c) 

— e = c = hf(e) + 7X0 - Xi + piz*, 


(28) 

(29) 


where h — 1, 7 I, pi = 1, = 1, and the coefficient of x\ is —1. The result 

is checked by the analog computer study shown in Fig. 4. 


Example 3: Given a control system shown in Fig. 1 with the linear part 
transfer function 



Fio. 4. Analog computer results — example 2. 


and /(e) = e + e* which is continuous but not bounded. Note that h = 1 and 
the system has one pole at — 1 and one pair of imaginary poles. The derivative 
is f'(e) — 1 + 3e* and 1 + hf'(e ) = 2 + 3c s > 2, and /( 0) = 0. From condi- 
tion Eq. 6, k has to be greater than 1 + e J . For e arbitrarily large, k has to be 
chosen as infinity. Then setting q to zero in inequality 21 gives, for l/k — 0, 

R* a.(M - Re [1 - > 0. (31) 

This system is asymptotically stable according to Corollary Ila. Referring to 
Fig. 3, y *= 0, h == I, pi = 1, and wi = 1, but the (s + 1/3) part is replaced 
by + 1), with coefficient —1 for x \ . The system equations are: 


(32) 


it = — x t +f(e) 

i\ — 

i\ *■ — z* +/(e) 

— « = c — /(c) — ii + z*. (33) 

This result is checked by the analog computer study shown in Fig. 5. 

Example 4: Given a control system shown in fig. 1 with the linear part 
transfer function : 

«•<*> - firr - rh + 1 M 

and/(c) = c + e*. This differs from Example 3 in having a pole at —2 instead 
of at —1. Setting q = 0 and l/k = 0 in inequality 21 gives 

R.G.(W-K.[l-^]-|±^>0. (35) 



Fio. 5. Analog computer results — example 3. 


This system is asymptotically stable according to Corollary Ila. This result 
is checked by the analog computer study shown in Pig. 6. 


CONCLUSION 

This paper has given two new theorems for the stability of control systems 
with one nonlinearity /(e). These theorems are extensions of Popov's Theorems 
PI and PII. Theorem I removes the restriction that y > 0 so that y > 0 
and a pole at the origin may not be necessary. It also introduces a non-negative 
constant h in the expression of the linear part transfer function defined 

in Eq, 7. The inequality 8 in Theorem I includes a non-negative constant 
(l/k). Example 1 demonstrates its application to an asymptotically stable 
control system. Theorem II is applicable to control systems with its linear part 


transfer function G i (s) having pairs of imaginary poles. If the nonlinear func- 
tion f(e ) is bounded, it is shown that a similar test of inequality can be made 
for predicting the quasi-asymptotical stability of the system. Corollary Ila 
disregards the boundedness of f(e) and predicts asymptotic stability whenever 
a non-negative quantity q (used in Theorems I and II) can be set to zero. 



Fig. 6. Analog computer results — example 4. 


Theorem II is applied to Example 2 and Corollary Ila is applied to Examples 3 
and 4. The results cheek with those obtained from analog computer studies. 

APPENDIX I 


Lemmas 


(1) If G(s) is a real positive function of the form N(s)/D(s) where N ($) 
=* dn+i3 n + • • • + d 2 s + ch and D(s) = s n + • • * + a 2 s + <*i, then there exists 
a function «(t) with the following properties: 


f' **Mt > «,|y(D|* - 0i|y(e)l |y(DI 

(<)] = H (ju)&(jw) and I H(j») | 2 - Re [G(ja>)], 


where $ is the Fourier transform, oti and 0i are positive constants. (See (4), 
Lemma 1 and p. 15.) 


( 2 ) 


If p(0 satisfies the conditions Supo<t<r |p(OI < /*, J |p(0|d< < 


and J \dp(t)/dt\dt < p then p(l)<j>(t)dt 


< rip Supo<<<r|y(OI where p 


and ri are positive constants. (See (4), Lemma 3.) 

(3) If r 2 (t) < a Supo<«<r|r(<) | + 0, then |r(<)| < a + V/3, 0 < t < T 
where a and 0 are positive constants. (See (4), Lemma 2.) 


(4) If both J ef(e)dt and de/dt are bounded, then lim,-«e(0 = 0. (See 
(1), Appendix 5.) 

(5) If e(f) = 0 then limi-M, y (<) = 0. (See (1), Appendix 6.) 


Proof of Theorem I 


Let Gi (s) be expanded into the form 

Gi(s ) = h + y/a + (ci + c*« + • • • 4- c n s n ~ l )/D(s) (1-1) 

where D{s ) = a* + a„a n_1 + • • • + a*s + ax, a Hurwitz polynomial. This 
expansion is shown in Fig. 7. The system equations corresponding to Fig. 7 



m 


G,(sl 


<b) 


Fig. 7. Expanded block diagram of the linear part transfer function Oi(a). 


are given by Eqs. 4 and 5 where 


-►G 



r > 


ro 

0 • • 

• 0 1 


1 1 


r 7 ] 


x 0 


0 

0 1 • 

• 0 


0 


Cl 

X = 

Xi 

A = 

0 

0 0 I 

• 0 

b = 

0 

c = 

Cl 


X n 

W J 


0 

— ai • • • 

- a, J 


i 

♦ 

C n 


For the convenience of proving the new Theorem I, let the matrix A be simpli- 
fied to B by removing the first row and the first column. Similarly x, b, and C 
are simplified to get y, bi, and g. Thus, Eqs. 4 and 5 can be rewritten as 

y = By + bj/ (e) (1-2) 

*-/(«) (1-3) 

e = - rco - g'y - hf(e ) (1-4) 

where 


y = 

Xi 

Xi 

B = 

0 10- 
0 0 10 

0 0 0 1 

• 01 
• 0 
• 0 

b, = 

> 

•• © o 

J 

g = 

CO*’ 

v 


Xn 

w ■> 


— Oi 

-a„j 


w 


W 


and g' is the transpose of g. Let 



<t>r(t) = 0(t) =/[e(0] when 0 < t < T 

= 0 when t > T 

then the solution of Eq. (1-2) can be found in integral form as 

y(0 = W(0y(0) + J*W(t - u)bi<t> T (u)du, 


(1-5) 


where W(0 = £“‘C(sI — B) _1 ], the inverse Laplace transform is taken on 
every element of the square matrix. 

A new function jr(t) is defined as 


where 


jr(0 = — v(t) — qdv/dl — (h + 1/k + gy)0r(O 



g'W(/ — u)bi4> T (u)du. 


(1-6) 


(I-7> 


The Fourier transform of Eq. (1-6) exists and is found to be 

5[jr(03 — Jrijv) 


= ~ [(1 +i«?)g'(i<oI - B)-'b, + h + qy + 1 AXi«) 

= — [(1 + jo>q)Gi(ju) — y/jw + l/*>(j«) (I-8> 


where $(j«) is the Fourier transform of 4>rit). 

Let 

R(t) = J* j T {t)4>r(t)dt + J“ t*(t)dt (1-9) 

then it follows from Parseval’s theorem that 


R(t) = RelJrUu)**U<*)]d» 

= — J Re [1/fr + (1 + juq)Gi ( ju) ] | $ ( jw) | s dco 

+ ^ j Re [1/A- + (1 + ja>q)G i ( jcj) ] | <l> ( jw) | = 0. (1-10) 

This follows from the facts that l/k + (1 + sq)G t (s) is a real positive func- 
tion, by using Lemma (1) and Eq. (1-8). Substituting Eq. (1-5) into Eq. (1-4) 
gives 

e = - g'W(<)y(0) - v(t) - h<t> T (t) - yx 0 . (1-11) 


Differentiating both sides of Eq. (1-11) with respect to t gives 


1 Asterisk denotes the complex conjugate. 



« = - g'W(/)y(0) - t)(0 - h<t> T (t) - yx. (1-12) 

subtracting both sides of Eq. (1-11) by <j>T(t)/k, multiplying both sides of 
Eq. (1-12) by q and then adding them together yields 

— v(t) — qv(t) — (h + yq + l/k)<f> T (t) = [e — 4>r{t)/k 3 

+ foe) + g'[W(i) + ?W«)]y(0) + yx. + g**r(0. (1-13) 

It follows from Eqs. (1-6), (1-9), and (1-13) that 

R(T) = RAT) + RAT) + Ri(T) + RAT) + RAT) + RAT) (1-14) 

where 

Ri(T) = 1/k J\ke - «(0>(0di > 0 

since the integrand is always positive according to the assumption. 

/ r /•*(!*) 

m)di - q l (o) f(e)de 

“ q / /(«)* - q / JMde > - qF( 0) 

where «/• J* 


r*{*) 

F( 0) = / f(e)de > 0. 

R>(T) = f r g'[W(t) + qirmymm > - r ! I y (0) i Sup \ j(t)\ 

Jm * 0 <t<T 

by Lemma (2) since B has eigenvalues with only negative real parts. 

RAT) = J* yx<4{t)dt = 7 J\ oio* = 7/2[*o J (r) - **»«»] 

fl.(D = J T qh4>m{t)dl = gV2[^(D - * 2 (0)] > - 9^ 2 (0)/2 

MD = ^" «*(<)* > a,|y(r)|* - fii|y(0)| 17(701, by Lemma (1). 

Substituting all the R y s into Eq. (1-10) gives the following inequality: 
*i\j(T)\* + yx.*m/2 < |y(0)|CJi|y(DI + n Sup |y«)|] 

o<<< r 

+ 7.r o 2 (0)/2 + qF( 0) + qh^(0)/2. (1-15) 


r t = min. [on, 7/2] when 7 > 0 

= ai when 7 = 0 

— (81 + 7 i)/ r i 
r 4 * = l/2r, 


Now let 


then from the fact that | y (T*) < Supo<i<r|y(0 < Supo< ( <r|x(0l and 
Lemma (3), the inequality (1-15) becomes 

|x(<)| < r 3 j x(0) 1 + rt'Jyxo 2 (0) +2qF(0) + qhtfi 0). (1-16) 

This proves the stability of the system. 

Now from Eq. 5, we have 


e + hf(e) = - C'x. (1-17) 

Since x is bounded and /(e) has the same sign as e, both /(e) and e have to be 
bounded. Differentiating both sides of Eq. (1-17) with respect to t gives 

e = - C'x/[1 + hf (e) ]. (1-18) 

Hence e is also bounded because C'x is bounded and 1 + hf'{e) ^ 0. 

From the assumption of inequality 6, there exists a positive quantity 
N < k such that 


Ri(t) = l/k O - <H<)>(f)d« > [(* - A r )/A'] I c<t>(t)dt > 0. (1-19) 




/ T 

c4>{t)dt is bounded and this 

bounded value is independent of T. It follows from Lemmas (4) and (5) that 
lim^ e(i) = 0 and lim^ao y = 0. Again from Lq. (1-14) we have limi-*,, x 0 = 0. 
This proves the theorem. 



APPENDIX H 


Proof of Theorem n 


The transfer function Gi(s) is expanded as shown in Fig. 8. The system 
equations corresponding to Fig. 8 are given by 

i - Ax + b/(c), (II-l) 

Zu-i “ Zu (t “ 1# 2, • • •, m), (II-2) 

zu = — wfzu-i +/(c), (II-3) 

e = — C'x — p'z — hf(e), (11-4) 

where 


z = 

f ^ 

Z\ 

Z\ 

and p = 

0 ' 
Pi 


Z 7m 


0 



i 

. P m . 


and matrices A and b are the same as in Appendix I. 

Following the same procedure as in proving Theorem I, we have in this case 

R'{T) - R{T) + R 7 (T) + R*(T ), (II-5) 



f(e) 


— G 0 (s) 

“Irrar 




<w 


flO- 


G^s) 


(•) 


Fig. 8. Expanded block diagram of the linear part transfer function Gift) 
with poke along imaginary axis. 


where R(T) is the same as given in Appendix I and Ri{T) and R$(T) are 


RAT) = j p'z 4>(t)(H, 

(II-6) 

RAT) = q j T p'z 4 >(t)dt. 

(H-7) 





Using Eqs. (II-2) and (II-3), we can find Ri(T) as 
*t(D = E / PiZu<f>{t)dt = £ f PtZ z aZu + <*>» 2 z 2 «-ii 2 .-_i]df 

*— 1 Jo <— 1 Jo 

= 1/2 E P<[zu 2 (T) + c ofru-i'm - ZiH 0) - «M«-I*(0)D 


>a 2 |z(r)| 2 -0 s |z(O)| 2 , (II-8) 


where a* = min. [p;/2, p<w,*/2] and /3 2 = max. [p</2, p<w< 4 /2] for all t. Ob- 
viously, a 2 and /3 2 are positive constants. 

Assuming Supo<<<r !<#•(<) I = M, then R»(T) can be found as 


R*(T) 


= Q J' p'z = Q J <t>(0d( p'z) 

> -gM|p|[|z(T)| + |z(0)|] 


> - 2gM|p| Sup |z(<)|. (II-9) 

o<<<r 

Substituting Rt(T) and Rs(T) into Eq. (II-5) gives 

«il y(5^) i* + rzo s (T)/2 + as|z(r)| 2 < |y(0)|[j8,|y(T)| + r 2 Sup |y«)|] 

0<t£T 

+ 2g3f|p| Sup | z (t) | + 7Z 0 2 (0)/2 + qF( 0) 

o<«<r 

+ /i« 2 (0)/2 + ^|z(0)|*. (11-10) 


Now let 


-[:]• 


r 6 = min. [ai, 7/2, a 2 ] when 7 > 0 

= min. [a h <* 2 ] when 7 = 0 

r 6 = max. [20i /r 8 , 2r!/r 6 , 2|p|/r 6 ] 
r 7 l = l/2r* and r 8 = 20 2 . 

Then by the same reason as in obtaining inequality (1-16) we have 
l*l(7 7 )l* < r *[ | y (0) | + qM~\ Sup \ v (t)\ 

0<t<T 

+ r 7 2 [ 7 ^(0) + 2qF (0) + qh<f{ 0) + r 8 | z(0) |*]. (11-11) 

It follows from Lemma (3) that 
ll(0l < r«[|y(0)| + qAf] 

+ r,\W(0) + 2gF(0) + qh$ r ( 0) + r,|z(0)| 2 . (11-12) 

Inequality (11-12) shows that rj is bounded. Hence both x and z are bounded. 
In view of Eq. (II-4), e has to be bounded. Since 1 + hf'(e) j* 0, e is also 

bounded. For the same reason os in Appendix I, f c<t>(l)dl is bounded. It 
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follows from Lemmas (4) and (5) that lim^ e — 0 and y = 0. 

Now rewriting Eq. (II-4) as 

yx a + p'z = - e - hf(e) - g'y = € e (t). (11-13) 

Differentiations of Eq. (11-13) with respect to i can he written as 

(T^o) + (PlZt) + • * * + (p»Z2m) = €o (t) 

0 + Wl l (piZi) + + «*> m S (p w Z2m) = €l(<) ^ (II-14) 

0 + <*\ m {jp\Zt) + + Wm~~(pm*2m) = *m\t) 

The determinant of the simultaneous linear Eq. (11-14) is not equal to zero 
since w, j* «y ^ 0 for i j by assumption. The quantities at the right side 
of Eq. (11-14) approach zero as t approaches infinity. This shows that Eq. 14 
has a trivial solution as i approaches infinity. Thus, 

7Xo — PlZ 2 = P 2 Z 4 = * * * = PmZt m = 0, Or Xo = Z 2 = Z 4 — * ' * — 22m = 0 

as i approaches infinity. By Eqs. (II-2) and (II-3), Zi * z* = • • • = z 2m ~i = 0 
as t approaches infinity. 

This completes the proof of Theorem II. Corollaries Ha and lib follow 
quite obviously after the substitution of q = 0 into inequality (11-12). 
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